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7 )0 ۲ رو ۰( 
میلاد جهانگیری بهروز احمدی ندوشن 


چکیده تحفیفات زیادی در تکنیک شناسایی حسارت در سازه‌ها باستفاده از بهینه‌سازی تک‌هدفه انجام گردیده است. انحی را کاربرد 
روش‌های بهینه‌سازی چندهدفه موردتوجه محفقان قرا رگرفته است. از این رو در این پژوهش بااستفاده از الکوریتم‌های نوین بهینه‌سازی 
تکاملی چندهدفه و بیان توابع هلف سبختی و شکل مودی به‌منظور شناسایی خسارت در تیر دوسرساده و در تیر دوسرگیردار پرداعته شده 
است. نتایج حاصل از تیر دوسرگیردار با نتایج تجربی دیگر محقان ممایسه گردیده است. نتایج نشان می‌دهد که روش پيشنهادشده منجر به 


عملکرد مطلوب در موقعیت‌یابی و کیفیت‌سنجی خسارت می‌گردد. 


واژه های کلیدی شناسایی خسارت. بهینه‌سازی چندهدفه الگوریتم تکاملی. 
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معا صمتاهمنصستن۵ وتحصمتاص(۲0 ۳ 


محطوی0عآ- 01 طخ .ظ 1عصمطه[ 1۷۲۰ 


۵ ۱۳۱۱۵۱۱۲۵ ۵۲ ۱۵اه 0۵0 06وی هام۵۵ ۱۵۲ ۵9۵۵۱۵ ۸۵ 6اوها خر 
5 0۲ 60۱5۱6۲۵ 66۷ 0عاه6 ۲۱۵۷ موه 0۵000۲ 0۱۷۵ زرا0- ۱۱۵ 66۵ 1061111660110۲ 
۵ 01۱ 0190۲۱۱/۱5 0۵۱00 2۷۵1۵۵0 ۱۱۵-۵26۷۵ ۲۵۵۵ ۵ ۵0۴۵]25 0۲۱۵/6 11:5 .2۲۵۵36 
6 ۵7۱4 ۵6۵ ۱0۵0۲۱۵۵ اک 6 1 2611011 00۵8۵6 اک ۵ ک۵ 61۷0۵زر۵ ککت ای 0۵ 5/106 
6۵ ۱۷۲۱ 60۱۴۵۵۲۵۵ 0۳۵ ۳65۵/5 16 ,۳66 5۷/۴۵۴۵۵۳۱۵ 12060 16 ۳۵۲ ,۵۵ 1101۳16111 2601-6110 1۷۵۵ 
0 0 6101111066 8000 116 06111015۳۵1۵ ۳۵5۲۵/5 1116 .۲6۲۵/۱/۳۵ 1۱۱۵ ۲۵۲ کا5ع۳ 62۵6۲۱۱۱۵۵ 

0۲۵۵6 ۲ ۱06۵201101 ۳0۵ 0۵۵۱60101۷ 0 8۵ 


۱ 


* تاریخ دریافت مقاله ۹4/۱۰/۳۰ و تاریخ پذیرش آن ۹۵/۹/٩‏ می‌باشد. 
(۱) فارغ‌التحصیل کارشناسی ارشد مهندسی عمران دانشکده مهندسی عمران دانشگاه یزد. 
(۲) نویسنده مسئول: دانشیار مهندسی عمران, دانشکدة مهندسی عمران. دانشگاه یزد. 261 )2312011 002)طعه :تفر 
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مقد مه 
سازه‌ها در طول مدت زمان بهره‌برداری‌شان دچار آسیب 
می‌شوند. این عوامل به شکل‌های طبیعی و غیرطبیعی 
پومشازه انم کلارد غراف طتضی از قییل ولو بادای 
عوامل غیرطبیعی شامل انفجا ضربه و... تأثیر مخربی 
برروی سازه دارند. این اثرات مخرب گاهی برروی 
سطح خارجی سازه ظاهر می‌شوند که ترمیم و بازسازی 
آن به‌وسیلة بازرسی چشمی ممکن می‌باشد ولی گاهی 
این اثرات مخرب درون سازه رخ می‌دهد یا بنابه 
شرایط هندسی سازه قابل مشاهده به‌وسیله بازرسی 
چشمی نمی‌باشند. به‌علاوه. روش بازرسی چشمی برای 
ارزیابی شدت خسارت. ظرفیت محدودی دارد. بنابراین 
در چنین حالاتی تشخیص محل خسارت و همچنین 
شدت آسیب واردشده به سازه امری دشوار است و نیاز 
به انجام مطالعات گسترده و آزمایش‌های مهندسی دارد. 
تا افو ان رو ها سای وا اضعا 
در سازه می‌توان محل و شدت آسیب را شناسایی کرد. 
درنهایت می‌توان با بازسازی و ترمیم خسارت واردشده 
طول عمر مفید ساختمان را افزایش داد. به‌عنوان مثال 
می‌توان به حدود چهل درصد از زیرساخت‌های کشور 
آمریکا از قبیل پل‌ها؛ سدهاء راه‌آهن‌ها مراکز عمومی, 
بسا ستایهاه ادا راش وهای وت رم اصا ره موه که 
ازلحاظ سازه‌ای دچجار خسارت شده‌اند و نیاز به 
مقاوم‌سازی و ایمن‌سازی دارند [1] در سال ۱۹۷۵ 
برای اولین بار توسط هولاند تابع هدفی براساس 
فرکاشی وی رامق تام تیا ریت :هر اوه ات 
شد و بااستفاده از روش الگوریتم ژنتیک به‌صورت تابع 
تک‌هدفه به بهینه‌سازی آن پرداخته شد [2]. در سال‌های 
ل استا ومد ات کال سای سرا قدابحا یی تاد 
مطرح گردید که کاربرد آن به‌صورت تعریف معیار 
اطمینان مودال (۷۸6) مطرح شد [3] همچنین پاندی 
و بیسواز در سال ۱۹۹۶ معیار نرمی مودال را برای 
تشخیص آسیب در سازه مورد تحلیل قرار دادند [4]. 
سپس لیمیتره در سال ۱۹۹۱ معیار جدیدی برای 
۱[ 


انرژی کرنشی و معادلة مقدار ویزه بیان کول ۳ 
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سالاو در سال ۱۹۹۷ تابع تک‌هدفة تغییرات فرکانس 
طبیعی سازه را بررسی کرد [6]. معیار اطمینان مودال 
(0۷۸0 که توسط اوینز در سال ۱۹۸۶ اراشه شده بود 
برای یک مود ارتعاش از سازه کاربرد داشت و درنهایت 
گوا و اسپنسر در سال ۲۰۰۲ در همین راستا معیاری به 
نام معیار اطمینان مودال کلی (10۸) را ارائه کردند 
کشا مان مودهای تاکن مان فان فاد رذن 
درنهایت پررا و همکاران در سال ۲۰۱۱ معیار اطمینان 
مودال کلی را مورد بازبینی قرار دادند و معیار اطمینان 
مودال کلی اصلاح شده (۷11۷1201 را پيشنهاد کردند 


[7,8]. یی و لیو در سال ۰۲۰۰۱ درصدد بررسی عملکرد 
توابع چندهدفه براساس پارامترهای مودال برآمدند که 
با تعریف دو تابع تک‌هدفه شروع شد و به‌طورکلی 
تحقیقات خود را برروی توابع چندهدفه و کاربرد ایین 
توابع برای شناسایی آسیب در سازه‌ها شروع کردند [9]. 
در این پژوهش سعی بر آن شده است که با بیان تابع 
هدف جدید و استفاده از الگوریتم‌های بهینه‌سازی 
چندهدفةٌ نوین ازدحام ذرات و الگوریتم تکاملی مبتنی 
بر تجزیه به مطالعه و بررسی عیب‌یابی و تشخیص 
آسیب در سازه پرداخته شود [10-12]. درابتدا با 
مدل‌سازی یک تیر دوسرساده در سناریوهای خسارت 
یگانه و چندگانه بااستفاده از توابع هدف سختی مودی 
و شکل مودی و به‌کمک الگوریتم‌های بهینه‌سازی 
تکاملی چندهدفه به شناسایی و ارزیابی آسیب در این 
تیر پرداخته شده است. همچنین, نتایج حاصل از 
شناسایی و ارزیایی خسارت در مدل‌سازی تسر 
دوسرگیردار با نتایج حاصل از تحقیقات آزمایشگاهی 
هیو و همکاران برای بررسی پتانسیل توابع هدف و 
میزان عملکرد الگوریتم‌های تکاملی 10۳50 و 
( مقایسه گردیده است. 


فرمول‌بندی مسئلة تشخیص آسیب در سازه‌ها 
معادلة مقدار ویذه 
هد و ماک وکا انا ی ون 
آنالیز مودال سیستم‌های سازه‌ای نقش بسزایی دارده در 
موعالت»۱)سفادله تقدار ویو ساژه هن ال تا تالم و 


سال سیم شمارة یک» ۱۳۹۲ 


میلاد جهانگیری- بهروز احمدی ندوشن 


۲ معادلة مقدار ویِعءٌ سازه در حالت آسیب‌دیده 
موردتوجه قرار می‌گیرد. معادلة مقدار ویدهٌ سازه در 
حالت سالم و آسیب‌دیده به فرم زیر بیان می‌گردد: 


۰.۱۵ ,12,3 < 0:6 - ((۰۳(۰)۸((ه) - []) 
)۱( 
۱۸۵ .... ,1,2,3 ع 6 : 0 < (۰][(۰)۵8"(م) - []) 
)۲( 
که در اینجا ]۸۷ تعداد مودهای ارتعاش سازه 
است. ماتریس جرم در هر دو رابطه بدون اندیس 
می‌باشد. زیرا در هر دو معادله ماتریس جرم ثابت است. 
این فرض براساس مطالعات آزمایشگاهی می‌باشد. 
ماتریس جرم قبل و بعد از سیب بدون تغییر باقی 
می‌ماند و در این پژوهش تغییرات آن تابعی از خسارت 
افمالفله بهشازه تست در انتها ایس های ط 6 
به‌ ترتیب نمادهای سازه سالم. سازه آسیب‌دیده و شمارة 


شود ارتغافن ساره مرباشتد: [13] 


دسا ریت اوه 
احتمال وقوع انواع خسارت در سازه وجود دارد. گاهی 
اوقات به‌دلیل اجرای بد سازه از قبیل جوشکاری 
نامناسب. ترک‌خوردگی بتن ممکن است سازه دچار 
آسیب گردد. باتوجه به مطالعات انجام‌شده. ایین 
خسارت‌ها به‌طور مستقیم بر دو پارامتر سازه تأثیر 
می‌گذارد. خساراتی که باعث تغییر در مشخصات 
هندسی سازه می‌ گردد منجر به تغییر ممان اینرسی 
می‌شود و خساراتی که در فیزیک و مکانیک سازه تأثیر 
می‌گذارد باعث تغییر مدول الاستیک سازه می‌گردد. 
تغییر هرکدام از این دو پارامتر باعث تغییر در ماتریس 
سختی المان و درنهایت ماتریس سختی کل سازه 
می‌گردد. برای محاسبة ماتریس المان نیاز به تعریف 
شاحص خسارت المان می‌باشد. این شاخحص خسارت. 
پراساس مشخصات هندسی سازه (ممان اینرسی) و 
ات شک کی ‏ ق ‏ یو سک نی رس 
می‌ گردد. شاخحص مذکور میزان خسارت اعمال‌شده به 
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0 


محاسبه می‌گردد: 


0 
۳( (ظ۲۱ ,123 ع) ز لس - 21 رم 


که در اینجا تلط تعداد المان سازه می‌باشد. این 
شاخحص همواره بین ۰ و ۱ می‌باشد. صفر به‌معنای سالم 
بودن المان در سازه و یک به‌معنای زوال کامل المان در 
سازه می‌باشد [14]. 


سخنی المان آسیب‌دیده 
اکنون برای محاسبةٌ ماتریس سختی المان آسیب‌دیده 
پراساس علم مکانیک محیط های پیوسته خواهیم داشت: 


2 [رک1] ۰( ره > 1( 5 [1] 


در اين رابطه [,16۵] و [,16۵] به‌ترتیب ماتریس 
سختی المان 1 ام در حالت‌های سازه سالم و آسیب‌دیده 
می‌باشد [15,16] 


فرمول‌بندی مسئلة بهینه‌سازی چندهدفه 
معیار اطمینان مودال 
معیار اطمینان مودال برای اولین‌بار توسط آلمانگ و 
برون برای استفاده در آنالیز مودال مطرح گردید. این 
معیار توسط محققان بسیاری مورداستفاده قرارگرفته 
است. این معیار برای مشخص نمودن میزان هم‌بستگی 
بین دو بردار می‌باشد. معیار اطمینان مودال که در هر دو 
تابع هدف مورداستفاده قرارگرفته است در حقیقت یک 
معیار مبدل به بازةُ * و ۱ می‌باشد و برای نشان دادن 
رابطةٌ هر دو بردار دل‌خواه قابل کاربرد است (مقدار آن 
همواره بین ۰ و ۱ می‌باشد). در علم آمار می‌توان به نام 
ضریب هم‌بستگی از این معیار نام برد. در حالتی که 
مقدار آن برابر صفر باشد نشان‌دهندة عدم ارتباط بین دو 
بردار می‌باشد به‌عبارت‌دیگر هیچ گونه هم‌بستگی بین دو 


بردار وجود ندارد و در حالتی که مقدار آن یک باشد 


نشریه مهندسی عمران فردوسی 


11 


بدان معناست که دو بردار کاملاً به یکدیگر هم‌بستگی 
دارند و مستقل نیستند و بردار تغییر نکرده است. معیار 
برای دو بردار دل‌خواه ۷ , ۶ به شرح زیر تعریف 
کر داد [17]: 


/۷ 2 
۳۵۳۸[ (0۱ 


- ۱۸۸۷ 
معیار اطمینان مودال براساس سختی مودی 
تابع مدف اول بااستفاده از معیار اطمینان مودال و 
براساس سختی مودال تعریف شده است. ازآنجاکه 
سختی سازه بعد از خسارت کاهش می‌پابد و جرم هم 
ثابت می‌باشد بنابراین کاهش فرکانس طبیعی سازه یک 
امر بدیهی تلقی می‌گردد. به‌همین‌خاطر. نکتهة بسیار 
کلیدی در رابطه ماتریس سختی مودال در مورد ضریب 
فرکانس طبیعی می‌باشد که روند کاهشی آن مارا در 
کاهش ماتریس سختی مودال یاری می‌نماید و درزمينة 
بهینه‌سازی تابع هدف برای هم‌گرایی و سوق به جواب 
ضریب جرم مودال با تغییرات تابع هدف هم‌سو 
می‌باشد و مشکلات تغییرات کوچک در فرکانس طبیعی 
سازه را جبران می‌نماید. درحقیقت جرم مودال سازه بعد 
از آسیب در سازه افزايش می‌یابد و این افزایش منجر 
به کاهش بیشتر ماتریس سختی مودال می‌گردد و باعث 
خسارت قون کر۵۵ [8]. 


0( * [2 .|۳7 .(۲۰]۸ |۲5 - [1] 
2 ۸۱ (] ۱1۰ ۱۳۱۰ [1] 
123,۰0 6 :6۵07۰ :۵:7 :[۳] ث < ,[1] 
۸( 
0 ).۵0 ۷ ۳9( - 10] 


اکنون‌که ماتریس جرم مودال برای هر مود در دو 
حالت سازه سالم و سازة آسیب‌دیده محاسبه می‌ گردد. 


درایه‌های بردار جرم مودال همان درایه‌های قطر اصلی 


ماتریس جرم مودال در هر مود برای هر دو حالت سالم 
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و آسیب‌دیده می‌باشد. سختی به‌معنای میزان نیرو به‌ازای 
جابه‌جایی واحد می‌باشد و ماتریس سختی مودال هم 
به‌معنای سختی سازه تغییرشکلیافته تحت الگوی هر مد 
شکل می‌باشد. هر ستون از ماتریس سختی مودی 
به معنای ایجاد یک جابه‌جایی واحد در درجهة آزادی 
متناظر با همان ستون می‌باشد. اکنون درایه‌های این 
ستون میزان نیرو در تمامی درجات آزادی تحت شرایط 
ایجادشده می باشد. ازاین‌رو قطر اصلی ماتریس سختی 
مودال درحقیقت نسبت به خسارت یجاد شده بسبار 
گره جابه‌جایی واحد اعمال‌شده ر نان می‌دهد و 


درصورت وجود خسارت ما را بهتر آگاه می‌نماید. رابطة 


معیار اطمینان مودال براساس سختی به فرم زیر تبیین 


شده است: 


۱ 


: ههد )| ِ 
۱۳ 0 < ان طتممکگ) ۰ (تمم16)) 
ار 
معیار اطمینان مودال مذکور یک رابطه برای 
مشخص نمودن میزان هم‌بستگی بین دو بردار می‌باشد. 
این معیار هم‌بستگی همواره بین ۰و ۱ است که در 
حالت صفر به‌معنای عدم وجود هم‌بستگی بین دو بردار 
می‌باشد و در حالت یک به‌معنای تطابق کامل بین دو 
بردار است. بنابراین هرچه مقدار این عبارت به یک 
نزدیک تر باشد به‌معنای تطابق بیشتر دو بردار می‌باشد؛ 
انش کوب اسر بیس شا ری کنو 
شاخحص های خسارت به‌درستی تعیین شوند. همچنین 
4 از بای تک سنوی ات ازه قایتل کاززرد آدست 
درصورتی که در آنالیز دینامیکی براساس آیین‌نامه 
حداقل سه مود اول لازم می‌باشد. محققان برای برطرف 
نمودن این مشکل. معیار گفته‌شده را برای هر تعداد مود 
دل‌خواه محاسبه کرده‌اند. براساس علم بهینه‌سازی 
به‌جای ماکزیمم‌سازی از فرم استاندارد بسرای 
مینیمم‌سازی استفاده‌شده است. بنابراین؛ تابع هدف اول 
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میلاد جهانگیری- بهروز احمدی ندوشن 


يا همان معیار اطمینان مودال براساس سختی به شکل 
قر رنه رونت اشت: 
)۱۱ ۸0۵۲۴ ۳7[ - 1 < ,۳ 
معیار اطمینان مودال براساس شکل مودی 
تابع همدف دوم بااستفاده از معیار اطمینان مودال و 
براساس اشکال مودی برای شناسایی خسارت مطرح 
گردیده است. اشکال مودی سازه. معرف الگویی از 
تغییرشکل سازه در هر مود ارتعاش می‌باشد. پراساس 
مطالعات انجام شده مکان و محل خسارت در تغییرات 
بردارزهای انکال مودی تاش پسواین دازنشد 
به‌عبارت دیگر با اعمال خسارت در المان‌های مختلف 
بردارهای اشکال مودی الگوهای تغییرشکل سازه را 
متفاوت نسبت به قبل گزارش می‌دهند [18]. ازاین‌رو 
معیار اطمینان مودال براساس شکل مودی سازه برای 
عیب یابی و تشخیص آسیب در سازه گزينة مطلوبی 
می‌باشد. این معیار براساس رابطة زیر تشریح می‌گردد: 
,۱۷۸۵۳۲۲ 


: ۱ ز(معحصده)] ۵۵۱۰ )| ۳ 
رز( مودصدهظ)] ۰ (معهسیه ۱۳) ی (نطرعموط)] ۰( ۳۵)) 


)۱۲ 
معیار اطمینان مودال شکل مودی, رابطه‌ای برای 
مشخص نمودن میزان هم‌بستگی دو بردار شکل مودی 
سازهٌ سالم و آسیب‌دیده می‌باشد. اين معیار نرمالیزه‌کننده 
که همواره بین ۰و ۱ است در این تابع همدف 
درصورتی دارای مقدار صفر است که هیچ‌گونه 
هم‌بستگی بین بردارهای اشکال مودی عددی و 
آزمایشگاهی نباشد. بهعبارت‌دیگر مقدار صفر مبین 
تخریب و آسیب جدی در سازه می‌باشد و در صورتی 
این مقدار برابر یک می‌باشد که سازه کاملاً سالم باشد و 
هیچ گونه مر ندیده باشد. بردارهای عددی. بردارهای 
سازه در حالت سالم و پر از عدد می‌باشند درحالیکه 
بردارهای تجربی؛ بردارهای سازه در حالت آسیب‌دیده 
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۷ 


علت ماکزیمم‌سازی معیارهای اطمینان مودال سختی و 
شکل مودی بدان دلیل است که درحقیقت هم‌بستگی دو 
بردار را به بیشینه مقدار برسد. به‌عبارت‌دیگ هنگامی 
معیاره‌ای مذکور به‌عنوان توابع هدف اول و دوم 
به‌صورت هم‌زمان بهینه می‌گرداند که معیارهای اطمینان 
مودالشان ماکزیمم گردد و این درصورتی امکان‌پذیر 
است که شاخص‌های خسارت در المان‌های آسیب‌دیده 
به‌درستی محاسبه و پیش‌بینی گردد. علت اصلی این امر 
آن است که فرکانس‌های طبیعی سازه در هر دو حالت 
متفاوت می‌باشد و فرکانس‌های طبیعی آزمایشگاهی که 
کوچک‌تر از فرکانس‌های طبیعی عددی هستند هنگامی 
اجازه ماکزیمم‌سازی به معیار مذکور را می‌دهد که 
شاخص‌های خسارت به‌درستی تعیین گردد. از سویی 
دیگر تمامی محاسبات و بهینه‌سازی با عنایت به رعایت 
قیود مسئله و باز؛ تغییرات شاخص‌های خسارت و 
معیارهای مودال می‌باشد. مشکلات ارائه‌شده در تابع 
هدف اول از قبیل توانایی‌های کاربرد برای یک مود و 
غیراستاندارد بودن فرم بهینه‌سازی برای تابع هدف دوم 
هم وجود دارد. بنابراین رابطة تابع هدف دوم یا همان 
معیار اطمینان مودال براساس شکل مودی به فرم زیر 
ارائه ون کر وه 
۱۳۲۱ :۷۵۵۳۲۷۲ 7[ - 1 < و۳ 
متغیرهای طراحی و قیود مسئله 

به‌طورکلی فرمول‌بندی ریاضی مسئلهٌ شناسایی خسارت 
شرت ای کنو 

فز این تانظه داز کاعم ارت )»دشر 
خسارت هر المان را نشان می‌دهد. 


< [((,0۰))رمنسم۳] 
ب).. ) ((م)) 

۹13 
به‌علاوه مقادیر شاخحص خسارت هر المان 
محصورشده در باز ۰ و ۱ می‌باشد؛ لذاء قیود مسئلاة 

شتایایی ات هر این مفتوده بان ری کرفنت ط 


نشریه مهندسی عمران فردوسی 


#۸ 
(۱۵) (طالرت قرق1 کل 17 > 8 >:0 
۱۹7( رتیه ریه9 | < ررم۱9 


بنابراین» فلوجارت کین تشخیص آ کشت در سازه 


به شرح زیر می‌باشد: 


بیان مستله بهینه سازی چند هدفه 


۱ تعریف متغیر های طراحی (:.4) 


۲ تعریف توابع هدف ((:.۲,6۹ ول(زء)دظ) 


۳- تعریف قیود مسئله (1 عر.8 > 0) 


بهیته سازی چتد هدقه با استفاده از الگوریتم 
های تکملی ۳۲۳۴56 و انهن۷0۵ 


چند هدفه یا رعایت قیود یه 
جواب مطلوب رسیده است. 
به گونه ای که تابع خطا به 
مینیهم مقداز پرسد؟ 


شناسایی و ارزیابی 


خسارت در سازه 


شکل ۱ فلوچارت تشخیص آسیب در سازه 


مدل‌سازی عددی و آزمایشگاهی 
مدل‌سازی عددی تير دوسرساده 
در این پژوهش یک تیر دوسرساده تحت اعمال 
سناریوی خسارت تک‌هدفه و جندهدفه مدل‌سازی 
گردید. به‌طورکلی سه نوع سناریو برای اعمال به سازه 
وجود دارد. نوع اول آن سناریوی خسارت یگانه نام 
دارد و تنها یک المان دچار آسیب هی گرزهه نوع دوم [ 


سازه به‌شرح زیر می‌باشد: 
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شناسایی خسارت در سازه‌ها بااستفاده از الگوریتم‌های بهینه‌سازی... 


110 أ‌ 19 أ 18 أُ 17 ۱ 1 


۶۲1 | 22 | 3 ۳ 


07777773 ٩, 2۸ 


6 


شکل ۲ سناریوی خسارت یگانه 


شکل ۳ ستاریوی سارت چند گانه 


جدول ۱ مشخصات هندسی و خواص مکانیکی مدل 


طول تیر ۳ 7 
پهنای مقطع عرضی | 9 
ارتفاع مقطع عرضی | " ۳۰/۲ 

مسا مقطع عرضی ۸ ۰ 1۳۹ 
چگالی مصالح ۰ 18/9 
مدول الاستیک ۲ ۰۴ ۷۲ ۱ 
ممان اینرسی ۱ (0*۳۱/۱۲) 

طول المان ت "5 


02 ۳ 


0.25 0 


شکل ۶ مقطع عرضی مدل تیر دو سر ساده 


بهینه‌سازی جندهدفه بااستفاده از الگوریتم‌های 
تجزیه و با نرم‌افزار متلب انجام گردید. در الگوریتم‌های 
مذکور اندازه جمعست اولیه 0۰ درنظر گرفته شد. 
مدت‌زمان اجرای برنامةٌ بهینه‌سازی جندهدفه در یک 
سیستم پنتیمم ۱/۷۲ هرتز و با حافظهٌ ۱ گیگا بایت برای 
هر دو الگوریتم حدود ۱۱ هزار انیه به‌طول انجامید. 
نمودارهای بهينه جبهه پارتو برای هر دو سناریوی 
خسارت یگانه و چندگانه برای هر دو الگوریتم به‌شرح 
کش مر کرو 
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۵۲۱۲0۵0 اد۴0 ماعد۳ عام‌سنة ت 


12 « ۷۷۱ 


1 1 1 1 1 
4 265 26 261 26 2 27 271 22 
۲۱ - ۱۷۸۲۲ 10 


دا 


شکل ۵ نمودار بهینهٌ جبهه پارتو در سناریوی خسارت یگانه بااستفاده از الگوریتم 10۳50 


۵ ۶ ۲۳۵۵۸ ۳۵۳6۸۵ عاوسن۹ 3 


۷۷۱ 12 
۳ 
ء 
ط 
۱ 
1 


شکل 1 نمودار بهینهٌ جبهه پارتو در سناریوی خسارت یگانه بااستفاده از الگوریتم ۱/۱0۰ 


۶۵ ۲۳۵۵۸ ۳۸6/0 عامتات۱۲ 
1 ۲ 1 1 1 


۰ 00295 
ب ۰ : 1 : 1 1 0029 


یبآ روبق شیب بت یسح بت تم | سجن لس میب ابیت أوفوون 


هی و جک و مخت 3623 


12. ۷۷ 


ب ۰ ۰ ۰ ۰ ۱ 1 سرکجج00 


ک 
#۳ 1 1 1 1 ۲ 
0022 00215 0021 00205 002 00195 0019 020185 001 0.15( 
۴۲ - ۴1 


شکل ۷ نمودار بهینهٌ جبهه پارتو در سناریوی خسارت چند گانه بااستفاده از الگوریتم 0۳50 


سال سی‌ام؛ شماره یک. ۱۳۹۹ نشریة مهندسی عمران فردوسی 


۱ 
آ 


03 
۲1 - ۴ 


شناسایی خسارت در سازه‌ها بااستفاده از الگوریتم‌های بهینه‌سازی... 


1۲2 2- ۷۷ 


شکل ۸ نمودار بهینةٌ جبهه پارتو در سناریوی خسارت چندگانه بااستفاده از الگوریتم ۸/۲ظ1۸0 


ع طا عزدمصوه(1 2۲۱( 


۳۳ 
(1 0 


8 9 10 


5 6 
۴۱60۱6۵۶ 


(۸۱۱۸۵۵ ۲ 


1 2 3 1 


شکل ٩‏ عیب‌یابی و تشخیص آسیب در سناریوی خسارت یگانه باستفاده از الگوریتم 10۳50 


هرکدام از نقاط روی نمودار پارتو خود یک دسته 
توب برای شا عفن عای سارت المان‌های ساره را 
اقفر هل کی آز یهاش یی شاه بآ بنانتن 
وضعیت سازه همان تابع اختلاف یا تابع خطا در 
محاسبات می‌باشد. درحقیقت. در سناریوهای خسارت 
اعمال‌شده به سازه مدل‌سازی‌شده شاخحص‌های خسارت 
واقعی را داریم. شاخص‌های خسارت به‌دست‌آمده از 
روی نمودار بهينءة جبههپارتو را می‌بایست با 
ای های عبات قاقعی سایبه موی آنن فا یه 
به‌روش خطای مجموع مربعات بین این داده‌ها صورت 
می‌گیرد و از بین مقدارهای به‌دست‌آمده کوچک‌ترین 
جواب همان گزینة مطلوب یا جواب مناسب برای 
پیش‌بینی خسارت در سازه می‌باشد. رابطةٌ تابع اختلاف 
یا تابع خطا در محاسبات به شکل زیر فرمول‌بندی 


نشریه مهندسی عمران فردوسی 


گردیده است: 


642 ی > ,۳ ,(1 


محاسبات و تحلیل برای پیش‌بینی و تخمین 
خسارت در سازه مدل‌شده بااستفاده از هر دو الگوریتم 
صورت گرفته است و نتایج حاصل در نمودارهای زیر 
مشخصن می‌باشد, نکقة کلیدی که در تخمین آسیب: در 
سازه‌ها وجود دارد. ارتباط بین گره در المان‌های مجاور 
هم می‌باشد. المان‌های در نزدیکی و مجاورت المان 
سس وله دازام کره مشع ک هی باشتند و شاعص 
حسارت المان آسیب‌دیده هنگام تشکیل ماتریس سختی 
کل بر المان‌های مجاور تأثیر می گذارد و خطاهای 
کوچک اطراف این المان‌ها قابلاغماض می‌باشد. 
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میلاد جهانگیری- بهروز احمدی ندوشن 


۷۱ 


و ۱۱2 


(_ 09 


09 
كِ 0.8 
ابو 
كِ 6 ی 
۳ اوه 
اِ بو 2 
۳ 031 
كِ 02 
اس 01 
10 9 8 7 6 5 4 3 2 1 ۱ 
هه مها 
شکل ۱۰ عیب‌یابی و تشخیص آسیب در سناریوی خسارت یگانه بااستفاده از الگوریتم 05۸/۲ 
۱ 5او90و01۸ 1۳۲6۲ 
3 __1 / / / ۲ 
: ۳ 
لت ۰ 4 -01 
اب ب فا بِ 
۴ ۲۱6۳۵۴۵ 
شکل ۱۱ عیب‌یابی و تشخیص آسیب در سناریوی خسارت چند گانه بااستفاده از الگوریتم 10۳50 
1 واجمنوه01 1۳۲6۲ 
1 ۳ ] ۱ 1 
۳ 
۱ 5 8 ۱ 
۱ : متسد امه ۱ ۱ ۱ 
شکل ۱۲ عیب‌یابی و تشخیص آسیب در سناریوی خسارت چندگانه بااستفاده از الگوریتم ۱۸0۸/0 


مدل‌سازی آزمایشگاهی تیر دوسرگیردار 
معرفی مسدل تیسر دوسرگیردار._ در اینن پژوهش 
مدل‌سازی تجربی هیو و همکاران به‌صورت عددی در 
نرم‌افزار متلب مدل شد و نتایج حاصل از تشخیص 


سال سی‌ام شمارة یک» ۱۳۹۲ 


سب در سازه که به‌شکل یک سناریوی اعمال خسارت 
انجام‌شده است. در پایان گزارش گردید. این دسته از 
داده‌ها برای اعتبارسنجی پژوهش‌های تئوری و براساس 
نتایج حاصل از کارهای آزمایشگاهی دیگر محققان 
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۷۳ 


صورت می گیرد. هیو و همکاران در سال ۲۰۰۰ 
پژوهشی برای شناسایی خسارت برروی یک سازه انجام 
دادند. جزئیات مدل‌سازی به شرح زير می‌باشد: یک تیر 
دوسرگیردار به طول ۷۰ سانتی‌متر که به ۲۰ المان 
تقسیم‌بندی شده است. جنس مصالح (مواد) آلومینیوم 
می‌باشد که مدول الاستیک آن 002 2-70 و ضریب 
پواسون آن ۷20.3 و چگالی آن برابر 3ط/ع1 222700 


معی‌باشد. سسطح مقطسع مسدل اراهشده 


۲۱-00500060 می‌باشد. 


طول تیر 9 1 
پهنای مقطع عرضی | ۲ تس 
ارتفاع مقطع عرضی | ۱/۰۰ 
مساحت مقطع عرضی ‌ ۱۳ 

چگالی مصالح 0 "/ع۰ ۲۷۰ 

مدول الاستیک ۲ ترا ۱ ۷۰ 
همان ایترسی ۱ ۱/۱۲(۳۵) 
طول المان ع ۰/۰۳۳ 


۳ 


شکل ۱۳ تیر دوسر گیردار 


0.006 


0.05 ۰ 


شکل ۱۶ مقطع عرضی تیر دوسرگیردار از جنس آلومینیوم 


مدل‌سازی هیو و همکاران د رآزمایشگاه. تیر 
دوسرگیردار مذکور توسط هیو و همکاران در آزمایشگاه 
ساخته‌شلده است. خسارت اعمال‌شده به سازهٌ 
آزمایشگاهی به‌صورت فیزیکی و با اه اعمال‌شده است. 
در المان شمارة ٩‏ با اه سطح تیر از بالا و از پایین 


نشریه مهندسی عمران فردوسی 


شناسایی خسارت در سازه‌ها بااستفاده از الگوریتم‌های بهینه‌سازی... 


به‌اندازة یک‌چهارم مقطع بریده شده است. به‌عبارت‌دیگر 
نیمی از مقطع سالم مانده است و نیم دیگر مقطع از بین 
رفته است. همچنین مقادیر فرکانس‌های طبیعی سازةُ تیر 
دوسرگیردار در هر دو حالت سازهُ سالم و سازه 
آسیب‌دیده قبل و بعد از بریدن تیر با ارّه اندازه‌گیری 
قن انیت افذاژه کنر شرف کا تسه ان طیسی شا 
به‌وسیلة دستگاه اندازه‌گیری دقیق و سنسورهای 
محدودی برروی گره‌های المان‌های تیر صورت گرفته 
است. ازاین‌روی میزان اختلافات در فرکانس‌های طبیعی 
در هر دو حالت سالم و آسیب‌دیده در مدل‌سازی 
تئوری و تجربی امری طبیعی به‌نظر می‌رسد. تجهیزات 
دقیق اندازه‌گیری در آزمایشگاه هیو و همکاران و 
همچنین مدل شماتیکی از تیر دوسرگیردار آسیب‌دیده 
به‌وسیلة ارّه با مشخصات دقیق هندسی خسارت در 


تصاویر زیر ارائه گردیده است [19]: 


13322222 


0.225 0 


شکل ۱۵ تیر دوسر گیردار بریده‌شده با اه در آزمایشگاه توسط 


هیو و همکاران 


مدل‌سازی تشوری در نرم‌فزار. تیر دوس رگیردار 
ساخته‌شده در آزمایشگاه در نرم‌افزار متلب مدل گردیده 
است. مشخصات تیر مدل‌ساوي‌شنده دققا هنن 
مشخصات تير فوق‌الذکر می‌باشد. در این پژوهش, بر 
تیر دوسرگیردار مدل‌سازی شده در المان شمارء ٩‏ 
خسارتی به‌میزان ۰/۸۷۵ درصد اعمال گردیده است. 
له ای ۲ مق می افو هدزای ۲۶ الاو 
۱ گره است. در این سازه در هر گره دو درجه آزادی 
درنظر گرفته‌شده است و درنهایت ماتریس سختی سازه 
با اعمال خسارت در المان ٩‏ محاسبه گردیده است. 


سال سیم شمارة یک» ۱۳۹۲ 


میلاد جهانگیری- بهروز احمدی ندوشن 


1۳2۳1286 ۹۵۵0۵210-9 


5 | مس ِ ِ تا ِ َِ 


۳۵۲0۵۴ ۲201286-0.875 1عتآوو۸ 


شکل ۱۸ تیر دوسر گیرداز با اغمال ستازبوی محساوت 


جدول ۳ نتایج فرکانس‌های سازه برای تیر دوسرگیردار 


اع هس هه ده 


۷۳ 


تحلیل مودال سیستم. نتایج حاصل از تحقیق هیو و 
همکاران در قسمت سازه سالم بسیار مطلوب ارزیابی 
می‌گردد و نتایج بسیار به‌هم نزدیک می‌باشد. اگرچه 
نتایج دز مت سوه آسسب‌دیده کی اشتاات, واوزه انا 
قابل‌قبول می‌باشد زیرا مشاهدة چنین اختلافاتی به‌دلیل 
خطاهای آزمایشگاهی از قبیل عدم گیرداری کافی در 
تکیه گاه‌ها. عدم اعمال صحیح آسیب به سازه عدم 
استفاده از سنسور به تعداد کافی برای افزایش دقت در 
اندازه‌گیری فرکانس در هر دو حالت سالم و آسیب‌دیده 
و دیگر خطاهای موجود می‌باشد. 


نمودار بهینة جبهه پارتو. نکته کلیدی که در زمينة 
نمودار بهينة پارتو برای انتخاب گزينة مطلوب وجود 
دارد درمورد انتخاب یک جواب مناسب از بین یک 
دسته جواب بهینه می‌باشد. در این تحقیق از معیار حطا 
در محاسبات بسرای انتصاب گزينة مطلسوب از مان 
دسته‌جواب ارائه‌شده در نمودار بهینهُ جبهه‌پارتو استفاده 
شده است. 

در نموداری‌های بهینةُ جبهه‌پارتو به‌ترتیب نقاط با 
مختصات‌های ۰/۱۲۸ و ۰/۱۶۸ برای تیر دوسرگیردار 


بااستفاده از الگوریتم 10۳50 و ۰/۱۲۹ و ۰/۱۶۷ برای 


سال سیم شمارة یک» ۱۳۹۲ 


تحقیق عحاشنز هیو و همکاران [14] 
آزمایشگاهی عددی آزمایشگاهی عددی 
خسارت‌دیده سالم خسارت‌دیده سالم خسارت‌دیده سالم خسارت‌دیده سالم 
ت ت ش۱۵2۵ ۸/۳۳۲ ۸/۹۸۶ ۳۹+ ۱۷:۷۰ ۸5/۹۹ مود اول 
ِ ۳۳/۹۷/۵ ۱2 ۱۳۱۳۵۹۸/۹۲ ۱۱۳۹/۷۲ ۱۵ ۱۳/۹۷( مود دوم 
ِ بِ 2۳۱/۸۳۹۰ 2-۳۱/:۳ 2-۸" 2/۹۶ 2۳۱۳/۱۷۱۹ "2-۹۳ مود سوم 
تیر دوسرگیردار بااستفاده از الگوریتم 1608۸/۲ برای 


گر دیده‌اند. 


و ۱۲۴۳۹۵ ط< صعع۳ اهعهماط ۲۵۵ ۲۳۰-۲۵ که و0 ۲۳۵۶ م۲ اعهتاو0 


‌ 


12 ۷ ۸ 


۲۱ - ۴ 


شکل ۱۷ نمودار بهینة جبهه‌پارتو برای تیر دوسر گیردار بااستفاده از 
الگوریتم 0۳50 


مود ۱۲۵۲۸/۵ ۱۱۵ سم اعهملا ف«ط مت 1۳0 که محوهاظ توص واععوط و6 


۳ 


۴ 


شکل ۱۸ نمودار بهینة جبهه‌پارتو برای تیر دوسر گیردار بااستفاده از 
الگوریتم ۱۸0۴8۸/۲ 


نشریه مهندسی عمران فردوسی 


۷ 


شناسایی و ارژیابی خسارت در ساژه. بااستفاده از 
نمودار بهینة جبهه‌پارتو به بررسی خسارت در تیر 
دوسرگیردار پرداخته سشد. محور افقی نمودار تشخیص 
آسیب نشان‌دهندة شمارء المان می‌باشد و محور قائم 


شناسایی خسارت در سازه‌ها بااستفاده از الگوریتم‌های بهینه‌سازی... 


الگوریتم تکاملی چندهدفه مبتنی بر تجزیه بوده است؛ 
نتایج مطلوب‌تری را گزارش می‌کند. وجود خطاهای 
کوچک مجاور المان شمارة ٩‏ بدیهی می‌باشد. زیرا 
هنگام تشکیل ماتریس سختی کل سازه المان‌های ۸ و 
۰ از سناریوی سارت تأثیر می‌پذیرند و تغییر 
می‌کنند. نتایج حاصل از پژوهش هیو و همکاران با 
تایج حاصل در این تحقیق مقایسه گردیده است: 


موه ۷۲0۵۳60 16۵ صحءظ )م۱ ۲۱۵۱0 ۲۱۲۵۴۱60 ۵ )۲۵۱۵ مد ممتاه۱06۵)16] ءمهصدنرا تمس ت50 


1222 121 77 
۳۳۲۰۱۰۰ 
رصان :۱ ۰ ۱۰ 


20 19 18 و 16 15 14 13 12 11 


1 1 1 ۱ / / / 


۱۱ 
ه‌ 
1 


#۴ 9 


(جدمتای062)0186) )صصع۴۱ ۵۶ مات 


شکل ۱٩‏ مقايسة نمودارهای شناسایی و ارزیابی خسارت در سازهُ تير دوسر گیردار هیو و همکاران با روش الگوریتم جندهدفةً ازدحام ذرات 
1۷0۳۹0( 


)معا ۲۷۲6۵۲۸/۱ )۱ صفءظ1 )ما۷ 0 ۲۳۵-۳۵۵0 اه صم‌مماه] ل4صده صمنا 1۱0611 ععمصوظ امیت 8 
1 


۰۵ :271۰ 
(0۰ظ) ا 
(.صسن۵ ۲-1۳01۴۲۸۲ 


۱۱۱۱ 


شکل ۲۰ مقایسهٌ نمودارهای شناسایی و ارزیابی خسارت در سازهُ تیر دوسرگیردار هیو و همکاران با روش الگوریتم تکاملی چندهدفةٌ مبتنی 


بر تجزیه 0۴8۸/9 


نشریه مهندسی عمران فردوسی 


سال سیم شمارة یک» ۱۳۹۲ 


میلاد جهانگیری- بهروز احمدی ندوشن 


۷۵ 


عصصا) دموا ۱۷۲۱۵۲۸/۸ 0صد مقهرن ۱ مه صمعتهموصمت 
۱ ۱ ۱ ۱ ۱ 


:1 1 ۳ 1 1 
ت ۳ ا 
(1 .6۲ ۳/]۲0۳۴560] - 


11 12 13 14 15 16 17 18 19 20 


۱ 


10 


7 ااا را اک ۳( 


(۱۵۱۵۵۵۵۸) ۱۱۵۵ ۵۱۱۵۵۵( 
هو هه هم و 
4 زره طا ۳53 دیا 


ِ 
دا 


01 


1 ِ 


3 4 5 6 7 8 9 


(طمذای10691111) )صمصع۱ 0 انا ار 


شکل ۲۱ مقايسة نمودارهای شناسایی و ارزیابی خسارت در سازة تير دوسرگیردار باستفاده از الگوریتم چندهدفة ازدحام ذرات 10۳30 و 


الگوریتم تکاملی چندهدفه مبتنی بر تجزیه ۷10۸/۲ 


نتیجه گیری 

۱- شناسایی خحسارت در این پژوهش براساس دو 
رویکرد الف) تشخیص محل آسیب و ب) بررسی 
سل رخ آسیب دنبال گردیده است. شناسایی موقعیت 
هندسی خسارت و ارزیابی شدت خسارت در هر 
المان نیاز به دو تابع هدف قدرتمند برای عیب‌یابی 
دارد. ازاین‌رو بهینه‌سازی چندهدفه به‌منظور 
۹ ۳ آسیب در سازه گزینه مطلوبی تلقی 
می گردد, 

۲- تغییرات فرکانس طبیعی سازه تاثیر بسزایی در 
دارد؛ زیرا تغییرات فرکانس طبیعی به‌صورت جمع 
بنابراین» تابع همدف اول تابعی حساس نسبت به 
شدت خسارت در سازه ارزیابی می‌گردد. 

۳ مکان منیا در سازه در مقادیر بردارهای شکل 
مودی تأثیر بسیار زیادی دارد. به‌عبارت‌دیگ معیار 
اطمینان مودال براساس شکل مودی نسبت به محل 


سال سی‌ام شمار یک» ۱۳۹۲ 


آسیب‌دیدگی عکس‌العمل نشان می‌دهد. بنابراین؛ 
تابع همدف دوم تابعی مطلوب برای تشخیص 
موقعیت هندسی آسیب در سازه می‌باشد. 

6- مدت‌زمان اجرای الگوریتم بهینه‌سازی چندهدفة 
ازدحام ذرات برای ترسیم نمودار بهينة جبهه‌پارتو 
برای هر تکرار حدود ۰ تانبه به‌طول انجامیده است» 
درحالی که مدت‌زمان اجرای الگوریتم تکاملی 
چندهدفة مبتنی بر تجزیه برای ترسیم نمودار بهينهة 
انجامیده است. اگرجحه الگوریتم‌های بهینه‌سازی 
جندهدفه مورداستفاده در این پژوهش دارای نرخ 
تکاملی چندهدفه مبتنی بر تجزیه دارای عملکرد 
بهتری ازنظر سرعت و دقت نسبت به الگوریتم 
تکاملی جندهدفةٌ ازدحام ذرات دارد. 


نشریه مهندسی عمران فردوسی 


۷۹ شناسایی خسارت در سازه‌ها بااستفاده از الگوریتم‌های بهینه‌سازی... 
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۷۸ شناسایی خسارت در سازه‌ها بااستفاده از الگوریتم‌های بهینه‌سازی... 


تشه مهنداییی عم رآ فرخونی سال سی‌ام» شمارةٌ یک؛ ۱۳۹۲ 


